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Abstract 

G\ ' 

ON ' We study a Schrodinger particle in an infinite spherical well with an oscillating 

wall. Parametric resonances emerge when the oscillation frequency is equal 
to the energy difference between two eigenstates of the static cavity. Whereas 
an analytic calculation based on a twodevel system approximation reproduces 
£N| \ the numerical results at low driving amplitudes e, we observe a drastic change 

\ of behaviour when e > 0.1, when new resonance states appear bearing no 

j> ' apparent relation to the eigenstates of the static system. 
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We study in this article the behaviour of a Schrodinger particle confined in a spherical 
cavity with an oscillating boundary that constitutes a particular kind of time-dependent 
Q\ ■ perturbation. Our study provides a conceptually simple "laboratory" in which the subtle 
and nontrivial aspects of the resonant coupling between the oscillating wall and a particle 
trapped inside the cavity can be investigated. Our original motivation in this work comes 
from our attempt to construct a dynamical bag model of hadrons [[!J ; however, our results 
may bear implications on the physics of a wide range of systems such as cavity QED || and 
perhaps even sonoluminescence 0. 

The system of a one- dimensional vibrating perfect cavity with quantized electromagnetic 
fields has been well studied 0. It was found that the electromagnetic field energy density 
inside a cavity vibrating at one of its resonance frequencies concentrates into narrow peaks 
regardless of the detailed trajectories of the oscillating cavity wall Furthermore, the 

amplitudes of these energy wave packets grow rapidly in time, producing sharp and intense 
pulses of photons. The distortion of the vacuum fields arising from the cavity wall motions 
leads to dynamical modifications of the Casimir effects |7]], which represents a fundamentally 
important and interesting feature of quantum physics. The problem of a quantum particle 
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in a box with moving walls has also been studied with an analytical approach [10], but the 
possibility of resonances was not discussed, which is the main interest in this work. 

If the oscillation amplitude eRo is small compared to the original cavity radius Rq, 
perturbation theory can be used to calculate the transition amplitudes between two states 
of the unperturbed system. This corresponds to what is usually observed in experiments. 
However the non-perturbative solutions of the complete time-dependent Hamiltonian (H = 
H + Hi(t)), where H is the time- independent part of the Hamiltonian, can in principle be 
remarkably different from the perturbative ones and can give rise to non-trivial features. 

We consider, as a first step, an infinite spherical well with oscillating walls: 
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where R(t) = Rq(1 + esinz/t) = Ro/a(t). Transforming to a fixed spatial domain via 
y = a(t)r, y = \y\ < Ro, and renormalizing the wavefunction (f)(y,t) = ct~ 3 / 2 (r, t) in 
order to preserve unitarity, we have 
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can be considered a small time-dependent perturbation if e and v are small enough. 

Since Hi(t) commutes with L 2 and L, we can look for solutions that are eigenstates of 
the angular momentum. This allows us to separate the angular dependence from the radial 
one in Eq. to obtain: 
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Using first-order perturbation theory, one can easily calculate the coefficients of the 
solution's expansion in terms of the unperturbed eigenstates. If the initial state is chosen to 
be \i >= \n = k,l = > ((f) n ,o = V2nwj (niiy)), we have 
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cl(t) = l5 nk E k tidt'{l-a 2 (t'))- 
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The term due to ^-f^f is exactly canceled out by the diagonal contribution of —^^y-p. 
The last integral is analytically solvable for v = io nk = (E n — E k ) /% and yields 
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sin ui nk t — - ln(l + e sin ui nk t) 
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The secular term u nk t/e in Eq. |]is a typical sign of a resonance. Notice that the secular 
term does not multiply a periodic function and the amplitude e that we suppose to be 
small is at the denominator. We can easily check that this is not a problem if we make a 
Taylor expansion of arctan (e + tan (u nk t/2)) / \/\ — e 2 in powers of e near e = 0, since the 
zeroth-order term exactly cancels the secular term. However the increase of c l n (t) in time 
remains. 

We can now calculate the expectation value of any observable as a function of time. We 
define the following dimensionless quantities: 
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E = mR 2 E/h 2 , {1 
v = mR\v jU . 

The perturbative results are in excellent agreement with the numerical ones when the cavity 
is oscillating out of the resonances. For example, at v = 7, e = 0.01 the fluctuations 
of the energy (Fig. |T]) correspond almost exactly to those of 1/R 2 (t), as one can expect 
from a quasistatic approximation, even though our system is not quasistatic. Even at high 
frequencies such as at z> = 90, e = 0.01, the first-order perturbative results are still acceptable 
(Fig. |2]a). Notice that in this case the energy is shifted up slightly and its fluctuations in 
time are smaller. This is due to the fact that the system is no longer able to follow the 
fast oscillations of the walls, and consequently the fluctuations as well as the value of the 
r.m.s. radius R s =< (y/R ) 2 > 1//2 are suppressed slightly (see Fig. |]b). 

At resonances, the perturbative approach breaks down and gives only an indication that 
a resonance exists. In order to study these resonances we solved the Schrodinger equation 
numerically, using a unitary numerical algorithm ||. For v — E 2 — Ei, we calculated the 
expectational values of the energy U =< E > and R s , choosing |n = 1,1 > as the initial 
state. In^Fig. | we plotted the results for I = and I = 1 (Z = , E 2 - E x = 14.8044 ; I = 
1 , E 2 — Ei — 19.7444) and two different values of e. The values for v = 7 are also plotted 
for comparison. The drastic change of behaviour of the system at the resonant frequency is 
evident even for very small amplitudes such as e = 0.001. 

At resonances the maximum expectation value of the energy, U max , varies as a function of 
e because of the trivial adiabatic factor a 2 (t) and, more importantly, non-trivial excitation 
processes. In Fig. [|we show m&x[a~ 2 (t)U] vs. e. For very small e (e < 0.002), the perturba- 
tion is not strong enough and the probability of exciting the second eigenstate never reaches 
1. The expectation value of the energy saturates (and equals E 2 ) for 0.006 < e < 0.1. In this 
regime, the frequency dependence of the energy maxima is well fitted by a Breit-Wigner 
function: £/ max = E\ + C/\{y — uq) 2 + r 2 /4], and the width Y increases linearly with e up to 
e ~ 0.1. For e > 0.1, even higher states are excited. 

Projecting the numerical solution on the eigenstates of the static system we found the 
expected result that for e < 0.1, the resonant dynamics is dominated by the lowest two 
eigenfunctions. This fact allows us to study the resonating system as a two-level system. In 
this case the differential equations for the coefficients reduce to: 
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V n (t)c 1 + V 12 (t)e- i ^ t c 2 \ , (8) 



c 2 = - 



^i(t)e^ Cl + \/ 22 (t)c 2 



(9) 



where Vij(t) =< i\Hi(t)\j >. Using the fact that Ci(t) changes little in a period T = 27t/uj 2 i, 
we can average Eq. |8] and Eq. ^| over a period to cast them into two coupled first-order 
ODE's with constant coefficients HI: 



Y,W ijCj . 



(10) 



Neglecting higher order terms in e, we have 
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W n = W 22 = , (11) 
4 f 1 - VT^) 

w 2 i = -w 12 = ^ — J -uj2i = n . (12) 

The system can then be diagonalized easily, giving Ci(t) = cos Qt and C2(t) = sin fit. 

When e < 1 then Q ~ 2o>2ie/3 and the period of the resonance lim e ^ ^~r = 2tt/Q = 00. 
In the other limit when e — > 1 then f2 — > 4/3co>2i ; but in this case our assumption that q(£) 
changes little in a period is no longer true and the averaging method no more valid. In Fig. |5| 
we plot the expectation value of the energy U = a 2 (t)(Ei cos 2 fit + E 2 sin 2 fit) and compare 
it with the numerical results. For amplitudes 0.005 < e < 0.1 the agreement is excellent. 

The matrix Wij can be written as — ifla 2 , where 02 is the second Pauli matrix. It follows 
that the vector formed by the coefficients c\ and C2 behaves like the spinor of a spin 1/2 
particle in a magnetic field along the j axis: 

i^ = n^ 2 |*> ■ (is) 

Therefore, if the initial state of the particle inside the oscillating cavity is one of the two 
eigenstates involved in the resonance, which corresponds to an eigenstate of S z , the evolution 
of the system will be a precession of < S > around the j axis. On the other hand, if 
the initial state corresponds to an eigenstate of S y we will obtain a stationary solution: 
\\f/(t) >= e Tint \^(0) >, which translates to 



Mv,t) = J^P-e^ ittt e-^Vi(y) iie-'^Vafo) 



(14) 
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The wavefunction in Eq. |TJ is periodical with period T = 2it/uj 2 i: 

<P±(y,t + T) = e w <f ) ±(y,t) , (15) 



where 9 = -2ir[E 2 /(hcu 21 ) ± 4(1 - v / T 3 ^ 2 )/3e]. 

We calculated numerically the solution choosing as initial function one of the two of 
Eq. |I4] at t — 0, and in Fig. [5] we show the resulting U. Although a 2 (t)U(t) is not strictly 
constant its variation is considerably smaller compared to other solutions. It is remarkable 
that such a highly dynamical system can show a quasi-stationary behaviour. 

For e > 0.1 the two-level approximation starts to break down. For e = 0.15 the third 
and fourth eigenstates become as important as the first two, and even more states are 
involved as one increases e further. The behaviour of the system changes drastically for 
e > 0.1, and we even observe the emergence of several new resonances that seem to have no 
straightforward explanation in terms of the unperturbed eigenstates. In Fig. |6| we show the 
maxima of a 2 (t)U(t) computed numerically for several driving frequencies choosing as initial 
state \n — 1, 1 — >. The resonance at v — 0021 is indicated, and it is much broader and 
smaller in amplitude compared to the new non-trivial resonances. It is interesting to note 
that even at these new resonances, the coefficients of the expansion in the static eigenstates 
are still approximately periodic. It may be possible to understand these new resonances for 
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e > 1 by including a few more levels in the two-level approximation. However the complexity 
of the system in this case warrants further study. 

For e < 0.005 the two-level approximation fails again; it continues to give the maximum 
of the expected energy as E 2 , typical of two-level systems, while in the complete system the 
energy maximum decreases as e is reduced. Also, the two-level approximation gives a period 
of the resonance T r greater than that of the complete system. 

We emphasize that the resonances we studied here are caused exclusively by the motion 
of the cavity wall, since the system has no interaction with electromagnetic fields. Another 
interesting feature of our system is the independence of its dynamics on R except for the 
rescaling of the oscillating frequency. 

It is also possible to consider a real system, hence with the electromagnetic interaction, 
in which an "oscillating-cavity" resonance occurs but the Rabi resonances do not. In fact, 
to observe Rabi resonances we need a cavity with radius Rq such that the fundamental 
frequency of the electromagnetic field u = 2tcc/Ro is equal to the difference between two 
energy levels, E n — oc h 2 n 2 j2mR\. It is hence not difficult to choose an R such that 
the Rabi resonances are not excited. In practice though, maintaining a stable mechanical 
oscillation with frequencies higher than some MHz is difficult. 

For simplicity we have only considered a spherically symmetric cavity with perfect wall. 
However, we conjecture that the resonances should not be too sensitive on the symmetry of 
the perturbation and on the detailed shape of the potential as long as the matrix element 
V12 (see Eq. |9|) is different from zero. One possibility is to use a microcrystal of conducting 
material with separations between the levels inside the conduction band of the order of 10~ n 
eV (~ 100 kHz). Forcing the crystal to vibrate at one of the resonant frequencies should 
excite many of the Fermi level electrons, which decay by emitting radiowaves. A second 
way could be to use a system with several, almost equispaced, energy levels. At a resonant 
frequency the particle, an electron or a trapped atom for example, absorbs energy from the 
driving oscillation to jump from one level to the next one and so on, as long as the resonance 
condition v ~ E n+ \ — E n is satisfied. In this way the frequency of the emitted quanta can be 
higher than the oscillation frequency, making them distinguishable from the electromagnetic 
noise due to dipole radiation at the driving frequency. 

In a further study we will consider a system with many equispaced energy levels and 
analyze the increase in energy with time. Ideally from such a system one can get quanta of 
frequency much higher than the driving frequency, and this is a major difference compared 
to the cavity QED situation, where at resonances typically a great increase in the number 
of photons with the same frequency as the driving force is expected. 

We thank Dr. C. K. Law for his suggestion of the two-level approximation. This work is 
partially supported by the Hong Kong Research Grants Council grant CUHK 312/96P and 
a Chinese University Direct Grant (Project ID: 2060093). 
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FIG. 1. Expectation value of the energy U vs. time (scaled by the oscillation period) for v = 7 
and e = 0.01. The initial state is chosen to be In = 1, 1 = >. 
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FIG. 2. (a) Comparison of energy calculated with the perturbative and numerical methods at 
a high oscillation frequency {p = 90, e = 0.01); (b) r.m.s. radius < (y/Ro) 2 > 1 / 2 for high (y = 90, 
solid line) and low {y = 7, dashed line) frequencies, both with e = 0.01. 
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FIG. 3. Same as Fig. 1, but for (a) I = 0, £ = 14.8044 and (b) I = 1, i? = 19.7444 with 
e = 0.001 (dashed lines) and e = 0.01 (solid lines). The arrows indicate the two levels E\ and 
For comparison we also plotted the dependence for is = 7, which show up as flat lines near E\. 



P 



100 



80 



60 



40 



20 



o /=0 
□ /=1 
x 1=2 



□ □ □ □ 



□ □ □ 



o o o o o o o 



0.00 



0.05 



0.10 

e 



0.15 



0.20 



FIG. 4. Maximum expectation value of the energy U max scaled by the "trivial factor" 
a 2 (t) = (Ro/R(t)) 2 vs. driving amplitude e for v = E2 — E\. 
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FIG. 5. Same as Fig. 1, but calculated with the two-level approximation (solid line), for 1 = 0, 
e = 0.02, v = E2 — E\ = 14.8044. The numerical results are also shown as diamonds. The dotted 
line shows the numerical result for the special case when \i >= l/y/2(\l > +i\2 >) 
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FIG. 6. Maximum expectation value of the energy U max scaled by the "trivial factor" 
a 2 (t) = (Ro/R(t)) 2 vs. the driving frequency for I = 0, e = 0.1 (solid line), 0.15 (squared line), and 
0.2 (dashed line). Several nontrivial resonances appear. 
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